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Abstract 

The 1 So-pairing gap in semi-infinite nuclear matter is evaluated microscopical- 
ly using the effective pairing interaction recently found explicitly in the coordinate 
representation starting from the separable form of the Paris iViV-potential. In- 
stead of direct iterative solution of the gap equation, a new method proposed by 
V. A. Khodel, V. V. Khodel and J. W. Clark was used which simplifies the procedure 
significantly. The gap A obtained in our calculations exhibits a strong variation in 
the surface region with a pronounced maximum near the surface. 
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Recently a method was elaborated of explicit consideration of the gap 
equation for semi-infinite nuclear matter within an approach based on the concept 
of the effective pairing interaction. The use of a separable form of the iVW-potential 
simplifies the problem significantly making it solvable numerically in the mixed 
coordinate-momentum representation. In Ref. [Q] the effective pairing interaction 
in the 1 S'o-channel was obtained explicitly without any form of local approximation 
for the separable 3x3 form ||, |]] of the Paris potential f|. In this case the gap 
equation is reduced to a set of three integral equations in coordinate space with 
kernels which are expressed in terms of u- and f-functions obeying the Bogolyubov 
equations with nonlocal gap A. These equations have the integro-differential form 
that makes their solution rather complicate. In Ref. || a method was elaborated 
of solving such equations, but it proved out to be rather cumbersome. As it is 
well-known, a rather slow convergence is inherent to the usual iterative method of 
solution of the gap equation. Since each iteration requires to solve these integro- 
differential Bogolyubov equations, a huge cpu time is necessary to obtain a reliable 
solution in this way. 

To make the procedure simpler, a new method of solving the gap equation 
for the case of nonlocal interaction is used (we refer to it as KKC). This method 
was originally devised for infinite nuclear matter but it is valid for any situation 
provided the gap A is much less than the Fermi energy £p. The main computation 
problem of usual iterative scheme originates from a strongly nonlinear form of the 
integral gap equation with a lot of iterations necessary to obtain the solution. The 
main idea of the KKC method is that the nonlinear form is important only for the 
magnitude of the gap but not for its momentum dependence. As to the latter, it 
is determined by integrals over a wide momentum space which does not practically 
depend on A, provided the parameter A/e-p is small. The gap A was identically 
represented in [0 as a product A(k) = Apx(k) of the constant Ap = A(kp) and the 
"gap-shape" function x(k) normalized to x(^f) = 1- Then the initial nonlinear gap 
equation was changed by a set of two equations. The first one, for the gap-shape 
function x{^), is a linear integral equation which does not practically depend on the 
value of Ap and can be readily solved. If the function x(k) is known, the equation 
for the gap amplitude Ap is just an algebraic nonlinear one which can be solved 
straightforwardly with standard methods. 

In the zero approximation of the KKC method, the gap-shape function x {fy * s 
calculated with Ap = 0. Then the zero approximation for the gap magnitude Ap 
is readily found. At the next approximation, the value of Ap' 1 is used as an input 
for calculating x (&) an d the procedure can be repeated. In principle, the KKC 
method contains no additional approximations in comparison with the standard 
one. It just rearranges the iterative scheme, making the convergence much faster. 
Already the zero approximation of the KKC method gives, as a rule, sufficiently 
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accurate results [0, ||. 

The KKC method can be readily extended to the case of nonzero temperature 
T |, |. In this case the gap-shape function is practically T-independent. The 
reason is that at small temperature which is of interest (the critical temperature 
is T c ~ 0.5Ap) redistribution of particles occurs in very narrow region nearby the 
Fermi surface which does not practically influence x(k)- Especially simple result 
appears in the zero approximation of the KKC method: 

A(k,T) = A F (T) x (k). (1) 

In the semi-infinite matter an additional coordinate dependence appears, but 
all physical reasons of eq. (l)-type representation remain valid with the coordinate 
and momentum dependent, but T-independent, shape function x an d the only T- 
dependent factor Ap(T). The latter one can be found as the asymptotical value 
inside nuclear matter which coincides with that of infinite system and is obtained 
by a much more simple calculation. To find the shape function we go to the limit 
T —>■ T c where the gap equation becomes linear. Such a method is used in this 
paper. It turns out to be, in the case of semi-infinite nuclear matter, much more 
convenient than the direct solution of the gap equation. 



In symbolic notation, the gap equation has the form R [TO 



A(T) = V c p ff ^(T)A(T), (2) 

where Vg ff is the effective pairing interaction acting in the model space in which the 
superfluid two-particle propagator Aq is defined. 

In the zero KKC approximation the temperature dependence of the gap opera- 
tor A(k, T) in homogeneous nuclear matter can be separated from its momentum 
dependence in the form of eq. (1). The shape function x(k) of the gap operator 
can be considered as a T-independent one up to the critical temperature T c . If T is 
close to T c , the gap A is negligible, and the propagator Aq(T c ) in eq. (2) coincides 
with that of the nonsuperfluid system which will be denoted as A (T). Hence at 
T = T c eq. (2) is reduced to a homogeneous linear integral equation for x(k). Its 
eigenfunction corresponding to the eigenvalue A = 1 is the shape-function we search 
for. 

Let us consider now the gap operator for semi-infinite nuclear matter which is 
nonuniform in the x-direction. We assume that the temperature dependence can be 
separated from the shape factor in the form similar to eq. (1) : 

A(xi, x 2 , k\, T) = A F (T) x(aci, x 2 , k 2 ± ), (3) 

where k^ denotes the momentum in the s-direction which is perpendicular to the 
x-axis. We use here the effective interaction calculated in H for the separable form 
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[0, ffl of the Paris potential with form factors gi{k 2 ) (i=l,2,3) depending on the 
relative momenta of interacting particles. In notations of 0, it has the form: 

V^ s (x 1 ,x 2 ,x 3 ,x 4 ;kl,k'l) = J2 A ij( X i X ')9i(k 2 ± ,x)g j (k'l,x'), (4) 

ij 

where X=(xi+x 2 )/2, X'=(x3+x^)/2, x=x\—x 2) x'=xz—x^ and gi(kj_,x) stands for 
the inverse Fourier transform of the form factor gi{kj_+k 2 ) in the x-direction. It is 
obvious that the gap-shape factor can be also written as a sum 

X(xi,x 2 ;k 2 ± ) = ^2xi(X) gi(kl,x). (5) 

i 

After substitution of eqs. (3) - (5) into eq. (2) at T=T C we arrive at the following 
equation for the components x%- 

Xi(X) =J2j dX x dX 2 A a (X, X t ) Bf m (X 1 ,X 2 , T c ) Xm (X 2 ), (6) 



In i 

where 

B^ m (X u X 2 ,T) = J ^^ dx i dx 2 g!(k 2 ± ,x 1 ) g m (k 2 ± ,x 2 )x 

x A (X 1 +^,X 1 -^,X 2 +^,X 2 -^-kl,T). (7) 

As is known [IT|, for the two-particle propagator A (T) at T > the Matsubara 
technique leads to the same expression as at T=0: 

A (r u r 2 , r 3 , r 4 , E = 2fi; T) = ^ \ — Ny ^ A (n)^(r 3 ) y> A /(r 2 )<^,(r 4 ), (8) 

a,a' ^-£a-£a' 

with T-dependent occupation numbers 

^( T ) = i +e(e !-,)/T » ( 9 ) 

where \i is the chemical potential of the system. The summation in eq. (8) is carried 
out over the single-particle states A = {n,k±} with functions <p\(r) = y n (x) e tk±s 
and energies e\ = e n + kj_/2m. 

After simple transformations we finally obtain 



bi(x,x',t)=-y: I 



dk ± l-N ni {k ± ,T)-N n2 (k ± ,T) 

■ X 



(2tt) 2 e ni +e n2 +k 2 ± /m-2fi 

xgLnMiVCnM,*')' (io) 



where 

9 l ni nMi X ) = J gi(k 2 ± ,x)y ni (X+^)y n2 (X~)dx. (11) 

The effective interaction V^ s is defined in such a way that the model space 
involves only the single-particle states with negative energies. Therefore the sum- 
mation over ni,n 2 and integration over in eq. (10) is limited by the condition 

The above formulae determine the kernels of eq. (6) for the gap-shape function 
x(xi,x 2 ,k 2 L ). Solution of this equation should be substituted into eq. (3) together 
with the factor of Ap(T)=0 which can be found at Xi 2 —>- — oo. Therefore it coincides 
with that of infinite nuclear matter. 

Before going over to solution of the above equations for the case of semi-infinite 
nuclear matter, let us check the accuracy of eq. (1), i.e. the zero approximation of 
the KKC method, in infinite matter for the specific kind of the iViV-potential in the 
separable form f|]. It was used for studying pairing in neutron and nuclear matter 
within the Brueckner theory in Refs. |12|, |13[ . We present here some results of these 
calculations in the form which is convenient for analysis of the T-dependence of the 
gap-form function. For the separable potential under consideration the momentum 
dependence of A is given by coefficients Cj of the sum 

A(k,T)=J2C l (T)g*(k 2 ). (12) 

i 

In fact, it is sufficient to check the T-independence of the ratios of these coeffi- 
cients to each other. It results that the T-independence of the ratios of GijG\ and 
C3/C1 is really true for both neutron and nuclear matter for different densities up 
to T = T c within 2%-=-3%- accuracy. An example for nuclear matter at kp = 0.8 
fm _1 which corresponds to the maximum value of A is given in Table 1. 

Let us return to consideration of semi-infinite matter. Eq. (6) can be considered 
as a particular case of a more general set of homogeneous integral equations 

Xi (X) = A(T) ]T J dX'K^X^'-^XiiX), (13) 

j 

where Kij(X, X'; T) are the kernels of this equation for arbitrary T and A(T) stands 
for the eigenvalue. The critical temperature can be found from the condition that the 
minimum eigenvalue Ai(T c ) is equal to unit. In principle, such a way which involves 
solution of eq. (13) for a number of values of T is possible but rather cumbersome. 
In the situation under consideration, when the pairing in infinite nuclear matter 
exists, it is obvious that T c in semi-infinite matter is the same as in infinite one. 
Therefore it is much simpler to find T c for infinite matter and use it in eq. (6) (or 
eq. (13) at A=l, T=T C ) as an input. 
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Table 1 



Demonstration of the T-independence of the gap-shape 
function for nuclear matter at kp = 0.8 fm -1 . 



T, MeV 









0.355 


—0.0271 


0.25 


0.355 


-0.0271 


0.50 


0.355 


-0.0271 


0.75 


0.354 


-0.0271 


1.00 


0.353 


-0.0270 


1.25 


0.350 


-0.0270 


1.30 


0.349 


-0.0269 


1.35 


0.349 


-0.0269 


1.40 


0.348 


-0.0269 


1.45 


0.347 


-0.0269 


1.50 


0.347 


-0.0270 


1.55(~ T c ) 


0.348 


-0.0271 



Let us now describe briefly how to calculate the set of integral equations, eq.(6). 
The kernels K im (X, X') of these equations are given by folding of the coefficients 
An(X,Xi) of the effective interaction V^ s with those Bf m (Xi,X';T c ) of the prop- 
agator Aq (T—T c ). The first ones were calculated in || for the value of chemical 
potential fi = —8 MeV which simulates the situation in finite nuclei. In this case 
the critical temperature is T c =0.124 MeV. In fact, the interval {X} = (— 8fm, 8fm) 
was considered only because properties of are trivial outside this interval. So, at 
X, X' < —4 fm the magnitude of V^ s coincides practically with that calculated for 
infinite matter, whereas at X, X' > 4 fm it tends rapidly to the free T-matrix taken 
at negative two-particle energy E = 2/i. As it is shown in 0, at a fixed value of 
X Q =(X+X')/2 the effective interaction Vg ff (A, A') vanishes rapidly when the rela- 
tive distance t=X— X' is growing, so that integrals involving V^ S (X, X') can be cut 
at \t\ > 4 fm. That is why the interval (— 8fm, 8fm) for the effective interaction is 
sufficient for any calculations. As to eq.(6) for the gap-form function, a much wider 
A-space should be considered, {A} = (— L; n , L ex ) with a minimum value of L in ~ 40 
fm (the value of L ex = 8 fm is large enough). The reason is that at small value of 
A F < 1 MeV the correlation length £ ~ vp/Ap is very big (£ > 10 fm). The distance 
between the left cut-off L in and the point under consideration (say, A > — 10 fm) 
should exceed a correlation length for the effects of the left boundary not to distort 
the asymptotic behaviour of the gap-shape function inside nuclear matter essential- 
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ly. For X, X' < —8 fm we use the effective interaction V£g(t) calculated for infinite 
nuclear matter. It is obvious that it depends only on the difference t=X—X'. 

The propagators Bf m (T) are new ingredients of the problem. They can be readily 
calculated in accordance with formulae (10), (11) using the technique elaborated in 

§• 

Such a large X-space makes it difficult to solve eq. (6) in the coordinate space 
directly along the way which was used in Ref. |2| for a similar equation for V^ ff . It 
is more convenient to use the Fourier expansion of the gap-shape function in the 
interval (— L in ,L cx ): 

XiW = E xVn{X) (14) 

where f n (x) are sin(27m(X— X c )/L) and cos(27m(X— X c )/L), L=L in +L ex , X c = 
(L ex — L in )/2). The kernels Kij(X, X') are also expanded in the double Fourier 
series. Finally, we deal with a set of homogeneous linear equations for the coefficients 

X?: ' 

3 N 

x? = EE W'xf, (is) 

3=1 n '=l 

with the matrix K^J 1 ' having not very high dimension. In this calculation we used 
the interval of {X} = (— 40 fm, 10 fm) (X c = — 15 fm), the value of iV=101 being 
enough to obtain the accuracy better than 1%. When the coefficients x2 are found, 
the gap-shape function x(X, t, k\) can be readily obtained with use of eqs. (14) and 
(5). To obtain the complete gap at zero temperature A(X,t, kj_,T = 0), one has 
just to multiply x by a constant A F = 0.167 MeV found for infinite nuclear matter. 

To present the results in a more transparent way, it is convenient to calculate 
the zero moment of x{X, x, s) over the relative coordinates: 

Xo(X) = JdsJdt X (X + t/2,X-t/2,s) = J £xi(X). (16) 

i 

Here, in accordance with Ref. [Q, the normalization g(k 2 = 0) = 1 was used. Let 
us also consider the function 

XF (X) = E XiW m(k 2 = k 2 F (X)), (17) 



where k-p(X) = J2m(/i—U(X)) is the local Fermi momentum (kp(X) = for 
fj,—U(X) < 0). The corresponding moments of the gap are as follows: A (X) = 
A F Xo{X), A F (X) = A f xf(X). The latter function is very important for the in- 
terpretation of results. It gives approximately matrix elements of A for the single- 
particle states taken at the Fermi surface. Analysis of this approximation for the 
case of the semi- infinite geometry is given in Ref. [[J. It should be noted that this 
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type of the local approximation is commonly used for description of pairing in finite 
nuclei | l4fl . Its accuracy for neutron stars was examined in |Tf|]. It worth to stress 



that in our paper it is used only for a more transparent presentation of results which 
were obtained without this approximation. 

Similar calculations were repeated for the chemical potential /a = —16 MeV 
which is the "self-consistent" value for semi-infinite nuclear matter coinciding with 
that of infinite matter. In this case we have T c = 0.57 MeV and Ap = 0.94 MeV. 

Both functions, A (X) and Ap(X), are drawn in Fig.l for both values of /i. As it 
is seen, all the curves possess pronounced surface maxima. Such a surface behaviour 
appears due to a big surface value of the effective pairing interaction 0. 

For more convenient analysis of the //-dependence of the surface effect, the gap- 
shape functions xf{X) for both values of \x are drawn together in Fig. 2. As it is seen, 
the surface effect in A is more pronounced for /i = —8 MeV. Indeed, in this case the 
ratio of the maximum surface value to the asymptotic one inside the matter is ~1.8, 
whereas it is ~1.5 for /i = — 16 MeV. Besides that, the position of the maximum is 
closer to the surface (X = 0) for fi = —8 MeV than for fi = —16 MeV. 

It should be noted that the surface effect for A turns out to be less pronounced 
than the one for the effective interaction itself. This is a consequence of a strong 
surface-volume coupling which takes place in the gap-equation due to a big pairing 
correlation length which is significantly larger than the width of the surface layer. 
Such a coupling is inherent to a pure quantum calculation and can be partially lost 
when a local approximation is used. It suppresses partially the surface maximum in 
A and makes the //-dependence of the surface effect smoother. 

Such a dependence of the gap on /i originates from two reasons. The first one is 
a direct energy dependence of the effective interaction Vg ff taken at E = 2/x M . The 
second one is the strong momentum dependence of the form-factors g-i{k 2 ) taken 
at k 2 = kp(X) = 2m(fi—U(X)). Both effects cooperate in enhancing the surface 
values of the effective interaction and of A for smaller values of \i. However, there 
is an effect which works in the opposite direction making the /i-dependence of the 
surface maximum of A less pronounced. This is a decrease of Ap at smaller values of 
\fi | which results in growing the correlation length. The latter smooths the surface 
effect itself and its //-dependence. 

The energy dependence of pairing should be important for consideration of the 
nuclear drip-line where the chemical potential /i tends to zero. The analysis of 
situation at small values of /i is in progress. 

The authors thank V. A. Khodel for the advice to use the idea of temperature 
independence of the gap-form function to calculations for semi-infinite nuclear mat- 
ter. They wish also to acknowledge S. T. Belyaev, S. A. Fayans and P. Schuck for 
fruitful discussions. Two of us (E. S. and M. Z.) thank INFN and Catania University 
for hospitality during their stay in Catania. 
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Figure captions. 



Fig. I. The zero moment A (X) of the gap (dashed lines), and the function 
A F (X) (solid lines) for two values of /i. 

Fig. 2. The gap-shape Xf(-^) for //= — 8MeV (dashed lines), and for \x= — 16MeV 
(solid lines). 
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